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Measured magnetic moments of single-particle states in the Pb region differ significantly 
from their Schmidt values. We discuss the reasons for this in terms of meson-exchange 
"^q", currents, isobar currents, core polarisation and other effects. 

o\ 

^ ' 1 Meson-exchange currents 

Nuclear physics calculations are based on the nonrelativistic Schrodinger equa- 
C/} \ tion in which wavefunctions for a many-body state are computed as a first step and 

expectation values of observables calculated as a second step. For processes such as 
7- and /3-decay the observable is represented by one-body operators, whose particular 
form can easily be deduced once the relativistic interaction between the currents and 
the fields is written down. For example, the interaction of a nucleon charge current 
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CN ■ with an electromagnetic field is given by a Hamiltonian Ti. = —J^A^ where is the 

vector potential describing the field and J M is the charge current of a single nucleon: 

O . 
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Jpik) = iu{p) F i7m - ^M^uK 

where k = p' — p. Here u(p') and u(p) are plane-wave Dirac spinors, and F\ and F 2 
are the Dirac and anomalous coupling constants. A nonrelativistic form of the charge 
current J M = (p, J) is obtained by multiplying out the Dirac spinors with the Dirac 



X 



u(p), (1) 



matrices and keeping only terms to leading order in 1/M — this produces a one-body 
operator sandwiched between Pauli spinors: 
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p{k) = F! + C(l/M : 

= ^(p' + p) + ^^icTxk + 0(l/M% (2) 

where for economy the Pauli spinors associated with u(p') and u(p) have been omit- 
ted. A Fourier transform to coordinate space and a multipole decomposition!!! leads 
to the familiar one-body electromagnetic operators, whose matrix elements are then 
evaluated in the many-body system using, for example, shell-model wavefunctions. 
When the coupling constants for the nucleon current embedded in the nuclear medium 
are taken from the free-nucleon system then the procedure is called the impulse ap- 
proximation. It has been widely successful, but at some point it is known to break 
down. This is because the nucleons in the nucleus are interacting through the ex- 
change of mesons and the perturbing electromagnetic field can disturb this exchange 



and even interact with the exchanged meson itself. As has been dramatically shown 
by Kubodera, Delorme and Rhoi with soft-pion theorems, the most important meson- 
exchange corrections (MEC) occur in magnetic moments and transitions (operators 
originating in the space component of the vector current, J M ) and in axial-charge 
transitions in beta decay (operators originating in the time part of the axial current, 

There are different approaches to the construction of the two-body MEC opera- 
tors that are equivalent up to certain unitary transformations!. They are: (i) the 
quasipotential formalism in which one starts from the covariant four- dimensional 
Bethe-Salpeter equation and uses the Blankenbecler-Sugar reduction to get a three- 
dimensional equationH; (ii) the equations of motion method in which unwanted de- 
grees of freedom are eliminated by a unitary FST transformation!, and the non- 
relativistic reduction effected by a Foldy-Wouthuysen transformation!!!!'!; and (Hi) 
the extended S-matrix formalism in which the S-matrix is expanded in a series of 
Feynman diagrams0'0. Our calculations are based on this latter approach using a 
phenomenological chiral Lagrangian for the isovector mesons, tt, p, and and 
augmented by the standard isoscalar meson Lagrangians for a and u mesons. 



1.1 'External' corrections 

There are two types of Feynman diagrams to consider: (i) those in which the 
electromagnetic field interacts with the nucleon charge current, J^k), to be known 
generically as "pair" graphs; and (ii) those in which the electromagnetic field interacts 
with the exchanged meson's charge current, to be known as "current" graphs. There 
is one comment we wish to make concerning the construction of the "pair" graph. 
The general form of its expression involves the product 

J,(k)S(Q)T(-q), (3) 

where J M (fc) is the nucleon charge current, Eq. (p]), S(Q) is the nucleon propagator 
of momentum Q, and T(—q) is the meson-production vertex of momentum q. To 
avoid double counting, the contribution from the standard one-nucleon current must 
be carefully identified and separated from the "pair" graph leaving the genuine two- 
nucleon MEC operators. This separation is usually done by dividing the nucleon 
propagator into positive and negative frequency components 

S(Q) = S^(Q) + S^(Q) (4) 

and retaining only the negative frequency component, S^~\Q). This gives the correct 
result providing J^k) and T(—q) are energy independent, i.e. J^k) does not depend 
on ko and r(— q) does not depend on go- However, for the magnetic moment, which 
derives from the spatial component of J M (fc), Eq. fll]), there is a dependence on k 
coming from the anomalous term (F 2 /2M)o-ioko- In this case Adam et a/.0 define 



Uk) = J^ + koJ^k), (5) 
where J M (fc) has no dependence on k . Thus J'Jk) = d/ (dk^^J^k). Similarly defining 

S M(Q) = ±S'M(Q) (6) 

then an additional term from the "pair" graph that contributes a genuine MEC 
involves the combination 

Jl(k)S'( + \Q)T(-q). (7) 

Such corrections are called 'external' by Adam et a/.0. Similarly when T(—q) depends 
on go there are 'vertex' corrections. The external corrections contribute to leading 
order, 0(1/M 2 ), for scalar and vector meson exchanges, and to next-to-leading order 
in 0(1/ M 4 ) for pseudoscalar meson exchanges. Here M is the nucleon mass. Vertex 
corrections only contribute to next-to-leading order. We take this opportunity to add 
the 'external' corrections to our previously computedli! MEC corrections for magnetic 
moments in the Pb region. The correction does not have a large impact, being of 
order ~ ±0.08^ with a plus sign for a single-particle proton state with j = £ + i or 
a neutron with j = £ — |, and a minus sign in the other two cases. 

1.2 MEC corrections for states in Pb region 

The procedure for constructing the two-body MEC operators for a magnetic 
moment calculation is as followsEl: 

• Write down expressions for the "pair" and "current" Feynman diagrams in terms 
of Dirac spinors for the nucleons for the case when the Lorentz index on the 
charge current, J^k), is space-like; 

• Expand the Dirac spinors in powers of (1/M) and retain leading terms, to be 
denoted J(k); 

• Construct the magnetic moment operator \x = — |iVfc x J(k) \k^o] 

• Fourier transform to coordinate space to obtain the two-body operator, ^(r, R), 
where r = ri — T2 and R = ^(r! + r 2 ) are relative and centre of mass coordinates; 

• Compute matrix elements of this operator in a many-body system using, for 
example, nuclear wavefunctions from a shell-model calculation. 

We can simplify the last stage by only considering matrix elements of the magnetic 
moment operator in closed-shell-plus (or minus)-one configurations. For the impulse- 
approximation one-body operator the calculation reduces simply to the single-particle 



Table 1. Meson-exchange current corrections to magnetic moments in 209 Bi and 209 Pb. 
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matrix element for the valence orbital outside the closed shell (the Schmidt value). 
Recall that the impulse-approximation magnetic-moment operator is 

H = g L L + g s S (8) 

with g L = 1 and g s = 5.586 for a proton, and g L = and g s = —3.826 for a neutron. 
The magnetic-moment expectation value for a single particle of angular momentum 
a is then 




The reduced matrix elements are defined in the conventions of Brink and SatchlerE-a. 
For the two-body MEC operators, the calculation becomes one of computing two- 
body matrix elements between the valence nucleon and one of the nucleons in the 
closed-shell core and summing over all the nucleons in the core. It is useful to express 
the result of this computation in terms of an equivalent one-body operator 

A* cff = 9h,e«L + #S,effS + 9p,es[Y2, S] , (10) 

where g h<eB = g L + 5g L etc.. Here g L is the free-nucleon coupling constant of Eq. (||) 
and Sg h the calculated correction. 

In Table [L] we give some sample results for magnetic moments for a valence O/ig/2 
proton and 0ii3/2 neutron outside a 208 Pb closed-shell core. Harmonic oscillator radial 
wavef unctions, hu = 7.0 MeV, are used. To accommodate the role of short-range 



correlations in the nuclear medium, we have multiplied the two-body MEC operator 
by a correlation function g(r) chosen to be g(r) = 9{r — d) with d = 0.7 fm. We also 
include vertex form factors at all meson-nucleon vertices, taking care to ensure that 
the equation of continuity remains satisfied^. 

Let us look at the 6g L value. In the S-matrix approach, the pion-pair and pion- 
current diagrams give large contributions but there is a significant cancellation be- 
tween them. For example, for the h 9 / 2 proton at 209 Bi we compute 5g h = 0.303 — 
0.247 = 0.056 with vertex form factors and 5o L = 0.325-0.230 = 0.095 without. This 
latter value is very elose to the estimate first obtained by MiyazawB of 5g L = 0.1 
from pion-exchange currents. We find, however, that heavy mesons make a substan- 
tial contribution to 5g h . Our result in Table |l] is 5g L = 0.056 + 0.130 = 0.186 for 
a Qhg/2 proton and Sg h = —0.032 — 0.057 = —0.089 for a O213/2 neutron, where the 
first figure represents contributions from pions, the second from heavy mesons. As a 
consequence, the magnetic-moment correction from MEC is positive for protons and 
negative for neutrons. 

The first experimental indication that the proton g L value is enhanced in nuclei by 
about 10% over its free-nucleon value came in the measurement of Yamazaki et alx3 
of the magnetic moment of the 11~ isomer in 210 Po. Nagamiya and YamazakJl! later 
showed the enhancement to be a general phenomenon present throughout the whole 
mass region. A systematic analysis^ of all magnetic moment data in the Pb region 
produces as best-fit values: 5g L = 0.15 ±0.02 for a proton, and 5g L = —0.03 ±0.02 for 
a neutron. Our results in Table [I] are in reasonable accord with this expectation, but 
there are other corrections, core polarisation, isobar currents etc., to be considered 
before a serious comparison with experimental data can be made. 



1.3 Isobar currents 

The MEC currents constructed so far involve nucleons interacting with mesons. 
There is a further process to be considered, in which the meson prompts the nucleon 
to be raised to an excited state, the A-isobar resonance, which is then deexcited by 
the electromagnetic field. Such proceses as these are called 'Isobar currents'. In Table 
it is seen that, the isobar current mainly quenches the spin g s -f actor by about 10% 
over its free-nucleon value producing a magnetic-moment correction that is positive 
for protons in j = £±~ orbits and neutrons in j = t—\ orbits, and negative otherwise. 



2 Core polarisation 

The nuclear wave function in a nucleus such as 209 Bi is more complicated than a 
single particle coupled to a closed-shell core and will have components of 2p-lh and 
3p-2h configuration whose impact on the magnetic-moment expectation value can be 
estimated in perturbation theory. To first order in the residual interaction, V, the only 



contributions arise from particle-hole states coupled to the same angular momentum 
as the multipolarity of the operator, A; namely, A = 1 for magnetic moments. At 
the 208 Pb closed shell, there are two such particle-hole states: proton {h^^hg/^} and 
neutron (£[3/2*11/2)5 an d their contribution to the magnetic moment of 209 Bi to first- 
and second-order is given by the expression!!!: 

(a\\An\\a) = 2j2T a ^ + 2^T a ^— (11) 

where Greek letters represent the particle-hole coupled states, viz. \ a) = | (/i~ 1 p Q! )A) 
of multipolarity A and the following notation is introduced: 



T a = <0||M||(/£ 1 Pa)A) 

L a = -a- 1 ((/i~Va)A I V I (a- 1 a)A) 

A a!i = ((fr~y«)A I V I (hp%)\) 

B af3 = (0 I V I {h- l p a )\ {hfp p )X) 

e a = e ha - e Pa (12) 

with a = (2a + 1) 1//2 . Here a represents the valence nucleon outside the closed-shell 
core. The energy denominators, e a , are estimated from known spin-orbit splittings of 
—5.6 MeV for the proton /i-orbitals and —5.85 MeV for the neutron z-orbitals. The 
residual interaction, V, is taken as a one-boson-exchange potential involving n, p, u, 
a and A% mesons and, for use in finite nuclei, multiplied by a short-range correlation 
function as described by TowneJli. 

The first- and second-terms in V given in Eq. (|TT|) can easily be extended to higher 
orders in the random phase approximation (RPA). It must be stressed that these are 
not the only second-order terms; there are others that we will discuss further in Sec. |3|. 
However, these terms are the simplest to evaluate. The selection rule that limited the 
number of particle- hole states, a, to just two still applies to all orders in RPA. The 
generalization of Eq. (|ll|) is now 



(a|| A M ||a) = 2j2T a [I-(A- B)/e)$ L p /e p (13) 

a(3 



where the matrix [J — [A — B)/ e] is first constructed and then inverted. Here I is the 
unit matrix. 

In Table |] we give some results for the core-polarisation calculation, based on 
Eq. (p~3|) and denoted CP(RPA). The calculation based on Eq. (|TT|) leads to an al- 
ternating sign series: quenching in first order, enhancement in second order and so 
on. As a consequence, the RPA sum is between 30 and 40% less than the first-order 



Table 2. Contributions to the equivalent effective one-body operator from all sources in the OBEP 
model for a 0/i-proton and a Oi-neutron in the Pb region. 
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calculation alone. The calculation mainly impacts on the spin g s -factor reducing it 
by about 25% compared to the free-nucleon value. The magnetic-moment correction 
is positive for protons in j = i — | and neutrons in j = £ + ^ orbits, and negative 
otherwise. 

2. 1 Coupling to core-vibrational states 

For single-particle states in the Pb region, especially the high-spin states, their 
wavef unctions are not those of pure single-particle configurations, but have some 
admixtures with low-lying vibrational states in the neighbouring even-even nucleus. 
For example, the lowest 13/2 + state in 209 Bi is not simply a proton 2 13 / 2 state but 
has a sizeable admixture of h 9 / 2 X 3~. Hamamoto0 has estimated the impact of core 
excitations on the magnetic-moment expectation values in a particle- vibration model. 
The correction, denoted 'Vib' in Table |2|, is small, < 0.06/iat, in all cases except the 
proton z 13 / 2 orbital and, to a lesser extent, the neutron i 13 / 2 orbital. 

2.2 Relativistic correction 

In deriving the MEC operators, a nonrelativistic reduction is made of the meson- 
nucleon interaction. This reduction is in essence an expansion in terms of p/M = v/c, 
where p is the typical nucleon momentum and M its mass. Terms to order (p/M) 2 
were kept. The standard one-body operator, Eq. (|8|), has only been determined to 
order (p/M) and we will consider its extension to the next-to-leading order, 0(p/M) 3 . 
The magnetic moment operator then becomesE!! 



^H 1 -^)-^*-^*™}^ 1 -^) < i4) 

It remains to estimate the expectation value of (p 2 /2M 2 ) for which we use harmonic 
oscillator wavefunctions. Since a second derivative is involved, this expectation value 
is very sensitive to the choice of radial wavefunction. The main result, denoted by 
'Rel' in Table ^ is a small reduction in the calculated magnetic moment of less than 
5%. 

3 Higher orders 

There are other second-order core-polarization corrections not contained in the 
RPA series discussed in Sec. |2|. Typically these terms have an operator structure, 
VfiV, where the residual interaction first excites a single-particle state to a 2p-lh or 
3p-2h intermediate state where the one-body magnetic-moment operator acts, then 
a second interaction reduces the intermediate state back to a single-particle state. 
There are no selection rules to limit the number of intermediate states, so the calcu- 
lations are computationally time consuming and have onlv been completed without 
approximation in light nuclei. In heavy nuclei, ShimizuE3 has explored a closure 
approximation to estimate these terms. 

There is another set of terms of the same order, namely first order in V but 
fourth order in meson-nucleon couplings, that have an operator structure 1/M + M7, 
where M is the two-body meson-exchange current operator. Again there are no 
restrictions on the intermediate-state summations. The importance of these terms 
was first pointed out by Arima et a/.i!@, where they are known as the 'crossing 
terms'. Their contribution to the magnetic moment is of opposite sign to that from 
second-order core-polarization terms involving one-body operators and cancels a large 
part of it. 

In Table |2], we quote the results of Arima et a/.H'B, where the cancellation between 
the second-order core polarisation, CP (2nd), and the 'crossing' term, MEC-CP, is 
clearly evident, but in detail it varies from case to case. The quoted numbers are 
model dependent and parameter dependent. 

Finally, we bring together all the computed corrections from core polarisation, 
CP (RPA), meson-exchange currents, isobar excitations and together under 'Other', 
the coupling to core- vibrational states, the relativistic correction and the higher-order 
processes, CP (2nd) and MEC-CP. In Table § the results are expressed in terms of 
an equivalent effective one-body operator. Consider first the orbital Sg L value. Our 
results of 5g h = 0.106 for protons and 5g h = —0.064 for neutrons are not too far from 
the empirical values determined in a best-fit analysis of magnetic-moment data in the 
Pb region of Yamazaki et aZ.@of 0.15±0.02 and — 0.03±0.02 respectively. The spin 5g s 
values are dominated by core polarisation. Our results of Sg s /g s = —40% for protons 



Table 3. Sum of all corrections to magnetic moments and B(M1) in the Pb region (last line gives 
7r S'7/2 case & t 132 Sn) in the OBEP model compared with experiment. 
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and —42% for neutrons show a large quenching consistent with an empirical relation 
often used: g S cff = 0.Qg s . In Table |3| we compare these corrections with all known 
magnetic moment and Ml-transition data in closed-shell-plus (or minus)-one nuclei in 
the Pb region. We also list the calculated results obtained by Arima, Shimizu, Bentz 
and Hyuga (ASBH)ll from core-polarisation and meson-exchange current processes. 
The agreement between theory and experiment is very good and is within 0.15// n in 
all cases except one, the neutron i 13 / 2 state in 207 Pb, where the strong cancellation 
between core polarisation and MEC seems to be too severe. The ASBH calculation 
has less of a cancellation here. 

Lastly, in Table |3] we have included one result for the 132 Sn doubly-closed shell 
core. The magnetic moment of 133 Sb has recently been measured at the OSIRIS mass 
separator of the Uppsala University Neutron Physics Laboratory by Stone et al. 0. In 
terms of the calculated effective g-factors the results are very similar to 209 Bi. The 



computed magnetic moment is in good agreement with experiment. 
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